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Abstract 

We present the STAR preliminary results on mid-rapidity and low transverse momentum 
mean, standard deviation, skewness, and kurtosis of net-charge and net-proton distributions in 
Au+Au and Cu+Cu collisions at -\fsim = 200 GeV for various collision centralities. All the 
measured high moments of these distributions can be scaled by the number of participating nu- 
cleons, consistent with the soft process emissions. The ratios of fourth to second order cumulants 
of both the net-charge and net-proton distributions are consistent with models without QCD crit- 
ical point. 



1. Introduction 

One of the goals of the physics program at the Relativistic Heavy-Ion Collider (RHIC) is to 
locate a critical point in the QCD phase diagram yj . The QCD phase diagram is often plotted as 
temperature (T) vs. baryon chemical potential (//b)- Lattice Gauge Theory (LGT) calculations 
indicate that at small chemical potential (jt/B~0), corresponding to the cases for RHIC and LHC, 
the transition from hadronic to partonic matter is a cross-over [2]. On the other hand, model 
calculations suggest that when the //b is larger, the transition becomes first order in nature fH]. 
Therefore, one expects the existence of a critical point at the end of first order transition []]]. At 
RHIC we plan to scan |4|] the center of mass energy ( Vsaw) to vary yUB in order to search for the 
location of the QCD critical point. Several lattice calculations suggest the existence of the critical 
point for /jb > 160 MeV [5, 6], corresponding to a center of mass energy of ysAw < 20 GeV. 

One of the characteristic signatures of the critical point is an increase in fluctuations of var- 
ious event-by-event observables 10. Hi- The moments of the conserved quantities such as net- 
charge, net-baryon and net strangeness are related to respective thermodynamic quantity suscep- 
tibilities. The LGT calculations, indeed, have shown that the fourth moment of these distributions 
have large values or diverge at critical temperature near the QCD critical point II5L |9L 1 1 OL 1 1 ill . It 
has been proposed that higher moments of event-by-event net-charge and net-proton multiplici- 
ties are significantly more sensitive to existence of critical point compared to measures based on 
second moments 1 12]. The fourth moments of these multiplicity distributions are expected to be 



proportional to the seventh power of the correlation length 111 2H . It is expected that the evolution 
of fluctuations from the critical point to the freeze-out may lead to non-Gaussianity in the event- 
by-event multiplicity distributions. Thus, the kurtosis of multiplicity distributions (net-charge or 
net-proton) could be a sensitive observable in the search of QCD critical point. 

We present the results of the moments of net-charge and net-proton distributions at midrapid- 
ity from Au+Au and Cu+Cu collisions at y/s^N = 200 GeV in the STAR experiment at RHIC. 
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2. Analysis Details, Results and Discussions 

With the large acceptance and excellent particle identification capabilities, STAR provides the 
best environment for studying event by event fluctuations. The charged particles were detected in 
the Time Projection Chamber (TPC) and protons were identified by their specific energy loss in 
the TPC gas. Standard STAR track quality cuts were applied in this analysis. The charged particle 
tracks are chosen to be within the transverse momentum (pr) range of 0.15 - 1.0 GeV/c whereas 
for protons the pj range is taken to be 0.4-0.8 Ge V/c in order to have a clean proton identification 
and to reject secondary protons from interactions with detector material. The analysis is carried 
out for Au+Au and Cu+Cu collisions at y/s^N = 200 GeV for various collision centralities. 
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Figure 1: Left panel shows the dE/dx (specific ionization loss) vs. rigidity distribution of charged particles for 70-80% 
most central bin in Au+Au collisions at -Jst/N = 200 GeV. Net-charge (middle panel) and net-proton (right panel) dis- 
tributions for Au+Au collisions at ysjviv = 200 GeV for three different centralities. The charged particle measurements 
are for \t]\ < 0.5 and 0.15 < pj < 1.0 GeV/c wheres the protons are within |v| < 0.5 and 0.4 < pj < 0.8 GeV/c. 

The left panel of Fig.[T]shows a typical dE/dx vs. rigidity (defined as the momentum x the 
charge of the particle) distribution for charged particles as measured in the STAR TPC for 70- 
80% most central bin in Au+Au collisions at y/lmf = 200 GeV. The figure gives a clean proton 
identification for pj < 1 GeV/c. The net-charge (h + - h~) and net-proton (p - 75) distributions 
for Au+Au collisions at yjs^N = 200 GeV for three different centralities (0-5%, 20-30% and 
70-80%) are also shown in Fig. Q] For the charged particles the distributions are symmetric 
around the mean value whereas the distributions for net-protons seem to be skewed in peripheral 
collisions. These distributions are further analyzed in order to extract various moments. 

The characterization of a given data set may be obtained by analyzing the moments of their 
distribution. For a given distribution of N data points, the third and fourth moments are expressed 
in terms of skewness and kurtosis as given by: 

Skewness = . anc j Kurtosis = ——r - - 3, (1) 

cr- 5 cr 4 

where 6N = N - (N) and cr = sj((6N) 2 ) is the standard deviation. Skewness is a measure of the 
symmetry of a distribution whereas kurtosis measures the peakedness of the distribution. For a 
Gaussian distribution, the values of skewness and kurtosis are zero. A positive value of skewness 
means a long tail on the right side and a negative value means the left tail is dominant. A data set 
with positive value of kurtosis is peaked at the center whereas a negative kurtosis is obtained for 
a flatter top near the mean. 
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Figure 2: Mean, standard deviation, skewness and kurtosis of net-charge (top panels) and net-proton (bottom panels) 
distributions as a function of centrality for ysjviv = 200 GeV Au+Au and Cu+Cu collisions. HIJING results for Au+Au 
200 GeV collision are also shown. 



Fig.|2]shows the evolution of various moments of the net-charge and net-proton as a function 
of collision centralities (denoted by the number of participant nucleons, A/ part ) for Au+Au and 
Cu+Cu collisions at y/s^N = 200 GeV. Note that the distributions have not been corrected for 
reconstruction efficiency and acceptance. The error bars include both statistical and systematic 
errors obtained by varying track quality conditions. In going from peripheral to central collisions 
we observe that the mean and standard deviation increase smoothly for all cases. The skewness of 
net-charge distributions are close to zero and similar for all centralities, whereas for net-protons, 
the skewness values have a decreasing trend in going from peripheral to central collisions. For 
net-charge distributions, the kurtosis values are small and has a decreasing trend in going from 
peripheral to central collisions for Au+Au systems. For net-protons, the higher values of kurtosis 
for peripheral collisions compared to central collisions indicate the distributions are more peaked 
at the mean value for peripheral collisions [13]. The results from HIJING event generator for 
Au+Au at 200 GeV are shown by solid lines, which are qualitati vely similar to those of the data. 



Several theoretical calculations based on Lattice QCD |9, 10, 11] and other QCD models 111411 
have suggested that the ratio of quartic to quadratic fluctuations of the net quark number, defined 
as /?4_2 to be a probe of QCD critical point. It can be shown that: 

= <(&V) > _ 3 /«jv) 2 ) = Kurtosis x Variance. (2) 
' ((6N) 2 ) 

R4 2 eliminates the 1/N behaviour necessary for the central limit theorem and removes explicit 
dependence on centrality. The significance of R4 7. is that it is entirely dynamic and can be 
sensitive to the location of the critical point l9l. ll0l.lllL 

Fig.[3]shows the R42 values for net-charge and net-proton distributions. While the net proton 
values are close to unity, net-charge values are higher. The results are similar to those obtained 
from HIJING. While the HIJING event generator is not modelled on the existence of QCD critical 
point, its ability to reproduce the qualitative features of the data needs further investigation. 
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Figure 3: Z?4,2(= kurtosis X variance) of net-charge (left panel) and net-proton (right panel) distributions for ^Js^n 
200 GeV Au+Au and Cu+Cu collisions. Results from HIJING are shown as solid lines. 



In summary, we have presented a study of mean, standard deviation, skewness and kurtosis of 
event-by-event net-charge and net-proton distributions in y/s^N — 200 GeV Au+Au and Cu+Cu 
collisions from the STAR experiment at RHIC. The results are presented as a function of collision 
centrality for mid-rapidity and low pj ranges. In case of net-charge distributions, the skewness 
values are close to zero for all centralities indicating the distributions are rather symmetric around 
the mean and the R\ 2 distributions are consistent with the predictions from HIJING. For net- 
protons, the kurtosis distributions are observed to decrease with increasing collision centrality 
and the R4 2 is unity for all centrality bins. The values are in agreement with model calculations 
from HIJING which do not have the physics of QCD critical point. Detailed comparison with 
other event generators, such as UrQMD, is in progress 11511 . 

The present study at the currently available RHIC energies only probes the fig < 30 MeV 
region of the phase diagram, while most theoretical calculations expect the QCD critical point 
to exist around fig >160 MeV. The current study will help to understand the expectations from 
various physics processes to the forthcoming RHIC critical point search program [4]. 
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